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Gravity forces cause distortions in tubular gas lenses. A theory is derived 
here which yields excellent quantitative agreement with measured distortions 
for various tube lengths, diameters, and gases. It is shown that in a gas 
lens of optimum design the displacement of the optical center has a maximum 
at the end of the lens. The amount of displacement increases with the fourth 
power of the tube diameter and with the square of the gas pressure. 

I. INTRODUCTION 

If a cool gas is blown into a hot tube (Fig. 1), the gas heats up first 
at the Avail of the tube and remains cool longer at its center. The density 
therefore, is higher in the center of the tube and decreases toward the 
wall. The increase in density is accompanied by an increase in dielectric 
constant. In this way the gas acts as a positive lens. 1,2 

At the same time, however, the cooler gas tends to sink down because 
of gravity, thus causing an asymmetric density profile in a horizontal 
tube. 3 Though a simple approach already gives an estimate of this 
effect,' 1 a more rigorous theory is derived here using a perturbation cal- 
culation which determines the transverse convection currents from the 
unperturbed temperature profile and then uses the currents to correct 
the temperature profile. 

II. TRANSVERSE CONVECTION CURRENTS 

The tube walls are at a temperature T w and AT degrees warmer than 
the entering gas. Heat diffuses toward the axis and determines the tem- 
perature field. Using the coordinate system shown in Fig. 1, the tem- 
perature field may be approximated by 2 

T = TV - Ar[l - 2 X -±-^ + (^-^) a } _I/ ' ) (1) 
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Fig. 1 — Convection currents and temperature distribution in a gas lens. 

where a is the tube radius and s a decay length given by the formula 

a\„ (2) 



s = 



7.3a 



v,„ is the gas velocity along the axis and a the thermal diffusivity denned 
as the ratio of heat conductivity k to heat capacity: 



P('v 



(3) 



The heat capacity is written here as the product of density p and specific 
heat at constant pressure. 

The temperature is related to the density p and the pressure p by 
the gas equation 
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V = RpT. (4) 

The density determines the gravitational forces gp which drive the gas 
particles in the transverse direction. The transverse components of the 
velocity field v(.t, y, z) can be found from Newton's Law 

dv 
gp = grad v ~ "pV"v + p -^ , (5) 

where v is the kinematic viscosity determining the frictional forces. The 
acceleration is described by the total differential dv/dt and for the steady 
state takes the form 

dv 

jj = (V grad) v. (6) 

In the problem under consideration the gas may be treated as a 
quasi-incompressible (Boussinesq) fluid. That means that variations 
of density may be neglected, except insofar as they modify the action 
of gravity. Forming the curl of (5) therefore, yields 

curl (gp) = —vp curl V 2 v + p curl-rr- (7) 

Using (4) and rearranging (7) one finds 

. . /grad T grad p\ . , , . dv . 

g X V^-Tp -) = v curl curl curl v + curl "ST" (8) 

Here grad p/p can be neglected compared with grad T/T, and T in 
the denominator will be replaced by the mean (absolute) temperature? 
T e . Finally, by inserting (6) one finds 

-(gX grad T) = v curl curl curl v + curl (v grad) v. (9) 

To solve this equation, a tentative velocity distribution is introduced 
which represents the flow lines shown in Fig. 1. The unknown coefficients 
are chosen in such a way that the equation 

div v = (10) 

is fulfilled, which assumes that the gas is incompressible. Then the 
velocity components 
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Vx = -% (a 2 - x 2 - y 2 )(a 2 - x 2 - 5y 2 ) 

v u = -4 — A xy(a 2 - x 2 - if) (11) 

^io / 2 2 2\ 

i'z = — (a - x - y) 
a 

result which leaves only the coefficient v unknown, since the velocity 
v, is determined by the forced laminar flow in the tube. v is the vertical 
gas velocity at the tube center caused by the gravitational forces. It 
may be assumed to be much smaller than the longitudinal velocity v... 
Though v is a function of z the variation of v in the z-direction is neg- 
ligible compared to its variation in the cross-sectional plane and has 
to be considered only in the acceleration term where dvjdz occurs 
multiplied with the velocity v 0: . 

With these approximations, v can be determined by inserting (1) 
and (11) into (9) which yields 

ig% e -'-y = m°£v + 18»..%- (12) 

Third- and higher-order products of x and y are neglected in this equa- 
tion since they are only important at the wall of the tube and contribute 
little at the tube center. 

Equation (12) is a linear inhomogeneous differential equation in z 
with the solution 



a AT a 
v a = - 



v T n 48 s - q 
where 



u ), (13) 



(-—•-•- d4) 



32 



A discussion of (13) is postponed in order to proceed with the cal- 
culation of the lens disturbance by using the derived convection cur- 
rents to correct the temperature profile which, in turn, gives the density 
distribution and the lens profile. 

III. DISPLACEMENT OF THE OPTICAL CENTER 

The gravitational forces cause a continuous flow of cool gas toward 
the bottom of the pipe, which distorts the temperature profile more and 
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more in the way shown in Fig. 2. The growing temperature gradient at 
the bottom, however, will increase the heat diffusion toward the center 
and counteract the convection effect. The equality of both effects is ex- 
pressed by the equation 

aV 2 r = v grad T (15) 

which determines the actual temperature profile under the boundary 
condition that T = T w at the tube wall. 

Considering that the temperature function for axial direction is 
much less curved than the radial one, d 2 T/dz 2 may be neglected and 
(15) separated with respect to z.° This yields 



f d 2 T . d 2 T 

a[^~9 + -7-2 



1 2 

dx 



dT . 
dx 



dT . T - T w 

T~ + V t 

By s 



(16) 



where T — T w is an exponential function of z as already introduced by 
(1) for the undisturbed temperature profile. 

No straightforward solution of (16) is known. Assuming, however, 
that the gravity effect, to first order, tilts the temperature profile in the 
^-direction as shown in Fig. 2, the amount of this disturbance can be 
calculated. The assumption implies that by transforming T(x, y, z) 
into new coordinates 



$ = x - 8(T W - T); n = y; f =z 



(17) 



the undisturbed profile can be regained, which in the following is denoted 
by 0(£, v, f)- Since this is symmetric with respect to £, the corresponding 
transformation in (16) must generate a differential equation for 6 
which contains only even terms in £. The requirement that the odd terms 




<*(T w -0){ = 




Fig. 2 — The temperature function T(x) and its transformation into a sym- 
metric function 0(x — dT). 
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cancel yields the following equation for 8: 

2a5 Ui F-**+ 8 *-— T~ W (18) 

The locus of the minimum of the temperature T(x) is of particular 
interest, for this is the optical center of the distorted lens profile. Fig. 
2 shows that this center occurs at a distance 

d = 8(T W - 0) £ .o (19) 

below the tube axis. Using for the undisturbed temperature profile 
given in (1) and solving (18) for 8(T W - 0) at £ = yields 



d = 



v„ 



2"f + 8^ 
s a 



(20) 



By inserting (2) and (3) into (20) one finally finds 

d = 7fo£fr^- -'-">• (21) 

The diffusivity a and the viscosity v for perfect gases are related by 
Eukens formula 7 



- = 7(9 -s- 

v 4 \ c. 



(22) 



c, being the specific heat at constant volume. As Table I shows, the 
decay lengths s and q given by (2) and (14) differ very little. Since (21) 
is not defined for s = q it is more convenient to use the following ap- 
proximation for (21): 



d 750 a 2 T s 6 



(23) 



which is valid for z < 2sq \q — s\. 

In Fig. 3 the displacement of the center of the lens profile is plotted 







Table I 






Cp/Cr 


a/y 


q/s 


He 

N 2 

C0 2 

CH., 


1.66 
1.41 
1.31 
1.31 


1.55 
1.35 
1.30 
1.30 


1.03 
0.93 
0.89 
0.89 
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60 80 100 120 

DISTANCE FROM INPUT IN mm 



Fig. 3 — Displacement of the profile center in a tubular gas lens of %-inch i. d. 
for flow rates of 0.5, 1, and 3 liters per minute using air. (a = 0.22 env/s). Meas- 
ured data by DeGano. 8 



versus the distance from the tube input for flow rates of 0.5, 1, and 3 
liters per minute. A tube of f-inch diameter and 100 °C wall temperature 
is assumed. The gas enters at room temperature. The mean temperature 
during the process is assumed to be T = 50 °C. The gas is air with a 
diffusivity a = 0.27 cm 2 /s. 

All curves show a linear increase of the displacement at the tube input, 
determined by the transverse acceleration of the gas. Further from the 
input the displacement follows the exponential decay of the temperature 
profile. The maximum displacement occurs at z = s. Measurements at 
the end of a 16-cm gas lens using the mentioned parameters are in fair 
agreement with the theory. 8 

In Fig. 4 the displacement is shown for a tube of J-inch diameter 
and two different gases: C0 2 with a = 0.125 cm 2 /s and N 2 with a = 
0.25 cm 2 /s. The temperatures are the same as in Fig. 3. The flow rate 
is 1 liter per minute. In this case, data are available for various tube 
lengths. 3 They show an excellent agreement with the predicted behavior 
of d versus z. 

The focal length of the tubular gas lens has a minimum if the flow 
rate is chosen in such a way that s equals the tube length. The maximum 
displacement occurs at the end of such a lens and has the value 



qa 
2040 a 



1 ^ V* 



AT 

T ' 



(24) 



A more useful measure for the gravity effect is the distance D by 
which a light beam has to be displaced off the tube axis to pass the lens 
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without deflection. Integrating at x = D over the tube length L one 
finds D from the requirement that the total deflection cancels: 



I 



dx\ 



dz = 0. 



(25) 



The development of the (disturbed) temperature field T about the axis 
yields for small distortion 

[ [D - d(z)]e-' /a dz = 0; (26) 

Jo 



and finally, by using (23) one has 



D = 



2/, 



4 1- (— +1K 2 ' 
1 ga \s 



3000 



1 - e 



-/./« 



(27) 



In Fig. 5 the displacement D is plotted versus the flow rate for CO. 
in a 7-inch tube assuming the same temperatures as in Figs. 3 and 4. 
Data measured by Steier 3 show good agreement with the theory. For 
L > s 



D 



ga 



3000 



(28) 



is a good approximation. According to this formula, the optical center 
of a C0 2 lens of optimum design would occur outside the tube if the 
tube diameter is larger than 1 cm. 
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Fig. 4 — Displacement of the profile center in a tubular gas lens of %-inch i. d. 
for a flow rate of 1 liter per minute using 00= (a = 0.1 onr/s) or N- (a = 0.2 
cm*/s). Measured data by Steier. 3 
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Fig. 5 — Displacement of the optical center of a tubular gas lens using CO- 
and a tube of 7 inches long and Vl-inch diameter. Measured data by Steier. 3 

IV. CONCLUSIONS 

The calculations show that the temperature distribution in a gas- 
filled tube undergoes a distortion which increases with the fourth power 
of the tube radius. A square law dependence on pressure is predicted for 
the range of 0.05 to 50 atmospheres where the thermal conductivity is 
independent of the pressure and therefore, the diffusivity a « 1/p. 

As a measure of the distortion, the displacement of the effective 
optical center in a tubular gas lens is calculated. Using C0 2 at room tem- 
perature and a tube of 10-mm diameter at 100 °C wall temperature the 
optical center occurs at the bottom of the tube. 
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